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ABSTRACT 

Hydrodynamical instabilities are believed to power some of the small scale (0.1-10 pc) turbulence and 
chemical mixing in the interstellar medium. Identifying such instabilities has always been difficult but 
recent observations of a wavelike structure (the Ripples) in the Orion nebula have been interpreted 
as a signature of the Kelvin- Helmholtz instability (KHI), occurring at the interface between the 
Hll region and the molecular cloud. However, this has not been verified theoretically. In this letter, 
we investigate theoretically the stability of this interface using observational constraints for the local 
physical conditions. A linear analysis shows that the Hn /molecular cloud interface is indeed KH 
unstable for a certain range of magnetic field orientation. We find that the maximal growth-rates 
correspond to typical timescales of a few 10 4 years and instability wavelengths of 0.06 to 0.6 pc. We 
predict that after 2 x 10 5 years the KHI saturates and forms a turbulent layer of about 0.5 pc. The 
KHI can remain in linear phase over a maximum distance of 0.75 pc. These spatial and time scales 
are compatible with the Ripples representing the linear phase of the KHI. These results suggest that 
the KHI may be crucial to generate turbulence and to bring heavy elements injected by the winds 
of massive stars in H n regions to colder regions where planetary systems around low mass stars are 
being formed. This could apply to the transport of 26 Al injected by a massive star in an Hn region 
to the nascent solar-system. 

Subject headings: ISM: kinematics and dynamics, Magnetohydrodynamics (MHD), Instabilities, As- 
trochemistry 



1. INTRODUCTION 



O yer 50 yea r s ago , it was postulated bv lSpitze"n ()1954f ) 
and iFriemanl (|1954D that hydrodynamical instabilities 
may form in interface regions between the hot diffuse 
gas ionized by massive stars and cold dense molecu- 
lar clouds. In particular, these authors suggested that 
elephant trunk, or spike structures, which are widely 
observed in star-formi ng regions, may resu lt from the 
Ray leigh- Taylor (see Chandrasckhaj] I1961D instability 
(RTI). Another classical type of inter face instability 
is the Kelvin-Helmholtz instability (see [Chandra sekharl 
1196 fl ) which occurs in the presence of a velocity sheer 
across the interface and is characterized by a wave- 
like periodic structure. Both types of instabilities have 
been considered to play a dominant role in the inter- 
stellar medium, as a power source fo r small scale (0.1- 
10 pc) turbulence (lElmegreenl 12004ft and in the mix- 
ing of chemical elements ( Rov fc Kunthl IT991I ). Obser- 
vationally, it has been hard to confirm the existence 
of these instabilities. Although the observed sizes of 
eleph ant trunks ma tch first order theoretical models of 
RTI (lFriemanlll954[) . recent observations of the velocity 
fiel d in the Pillar s of C reation and the Horsehead nebula 
by iPoundl (fl99l and iPound. Reipurth fc Ballvl (p003l ) 
tend to discard the RTI hypothesis for the formation of 
these stru ctures (in s tead, the selective photodissociation 
model of iReipurthl (119831) is invoked ). More recently, 
iBerne. Marcelino fc Cernicharol (|2010l ) (BMC hereafter) 
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observed a periodic wavelike structure (the Ripples here- 
after), at the surface of the Orion cloud (Fig. 1) which 
appears to be compatible with a KHI. However, this work 
was mostly qualitative and lacked a detailed model to 
assess if the development of the KHI is possible in con- 
ditions as those found in Orion. In addition, this former 
study did not consider the possible effect o f magnetic 
fields -which are known to be strong in Orion (jAbel et al.1 
l200l - on the KHI. Finally, even if there is evidence 
for the existence of the KHI in the interstellar medium 
(ISM), it remains unclear over which timescale this in- 
stability may convert energy into turbulent motion of 
the gas and hence if it can actually play a role in mix- 
ing the hot and cold gas. On the theoretical side, ex- 
tensive numerical magnetohydrodynamical (MHD) mod- 
els for the KHI have been developed and successfully 
applied to ex plain observed phenomena in solar-system 
plasmas (e.g.. IMatsumoto fc Hoshinoll2004T ). In addition, 
IMatsumoto fc Sekil (|2010T ) have studied in great details 
the saturation of the KHI and its evolution into turbu- 
lence and mixing of the gas. In this letter, we perform a 
linear MHD analysis applied to the situation in Orion's 
Ripples, using physical parameters determined observa- 
tionally. We derive the key parameters that characterize 
the instability and use them to determine the timescales 
over which the gas becomes turbulent due to the satu- 
ration of the instability. We discuss these results in the 
context of chemical mixing in the interstellar medium 
and transport of 26 Al in the solar system. 

2. OBSERVATIONS 

Fig -Ushows the Spitzer Space Telescope (|Werner et al.l 
120041) mid-infrared (mid-IR) image of the Ripples in 
Orion, which have been attributed to a KHI by BMC. 
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Figure 1. Spitzer Infrared Array Camera (Fazio ct al. 2004) 8 
(im image of the Ripples in Orion which have been attributed to a 
KHI by BMC. A denotes the spatial wavelength of the structure, 
Li in is the distance over which the instability appears to be com- 
pletely linear. The orientation of V^, the phase velocity of the KH 
wave at the surface of the cloud, is depicted qualitatively. 



KHIs occur at the interface between two fluids flowing 
relative to each other. In the case of Orion, the two 
fluids are the H n gas and the neutral gas of the molecu- 
lar cloud, and the velocity sheer results from the cham- 
pagne flow created by the Hn region bursting through 
the parental molecular cloud. When growing, the KHI 
gives to the interface a wavelike structure which can be 
seen in Fig. [1] An important parameter of the KHI is 
its spatial wavelength A which is connected to the phys- 
ical conditions in which the instability occurs. Here, A 
can be measured directly from the image and is found 
to be A = 0.11 pc for a distance to the Orion nebula 
of 414 pc. This corresponds to a spatial wavenumbcr 
k = 2tt/A = 2 x 10 -15 m- x . 

3. LINEAR ANALYSIS OF THE KELVIN HELMHOLTZ 
INSTABILITY IN ORION 

3.1. Objectives and method 

Our goal here is to study from the theoretical point of 
view, and using realistic physical conditions, the stability 
of an H II / molecular cloud interface against the KHI. In 
particular we want to determine wether magnetic fields 
can play a stabilizing role. It is also of great interest 
to derive some of the key parameters, for instance the 
growth rate 7 and tte range of acceptable wavelengths. 
In order to do this, we perform a linear stability study 
which includes magnetic field, compressibility and an an- 
alytical velocity profile across the sheered layer (see be- 
low). This differs greatly from the pr eliminary study 
of BM C which relied on an ideal case t Chandra sekharl 
1961) of incompressible fluids with a discontinuous ve- 
locity profile and no magnetic field. The present study is 
performed in a two-dimensional slab geometry, described 
in Fig. [2] for the initial conditions. These initial con- 
ditions are maintained by the magnetohydrodynamical 
equilibrium. The density gradient and velocity gradi- 
ents are along the y axis. The velocity and density pro- 
files across the interface are of hyperbolic-tangent form 
(M iura fc Pr itchctt 1982). The velocity is oriented along 
the x axis. The magnetic field direction is inside the 
plane defined by y and z and its orientation is defined by 
the angle 8 between B and z. For this first analysis, we 
have not considered the azymuthal dependance for the 



orientation of B, because this would imply heavy com- 
plications in the solving of the MHD equations. For the 
adopted configuration, the MHD equations are linearized 
and a perturbed quantity / can be expanded as a plane 
wave in the form of f(x, y) = f(y) exp{i(k x x — tut)} . The 
linearized equations can then be solved, with boundary 
conditions in the y direction and a given wave number k x 
for the corresponding eigen value (angular frequency and 
growth rate) as an eigen value problem. This is described 
in mathematical terms in the following section. 

4. LINEAR MODEL 
4.1. Basic equations 
The basic MHD equations are 

dp 
dt 



£ = -V.(pV), 



(1) 



ay 1 R 2 1 

- = - ( V.V)V--V(P + _ ) + -(B.V)B,( 2) 



dP 

0* 



= -(v-v)P-rp(v-v) 



~dt 



-cV x E. 



with the frozen-in condition 

V 

E = x B. 



(3) 



(4) 



(5) 



where r = 5/3 is a polytropic constant. The mass den- 
sity p and the magnetic field B arc normalized by charac- 
teristic values of p and B , the velocity V by the velocity 
jump across the boundary, Vo, the pressure P by Bq/8it, 
the spatial scale by the initial shear width L, and the 
time by L/Vq. 

4.2. Linearization and solution 

We consider perturbed quantities from an equilibrium 
state as, 

V = V + 5v, 
P = Po + 5p, 
B = B + b. 

The perturbed quantities are expressed as a plane wave 
in the form, 5 A = A(jj) exp [i(k x x — u>t)], where A, k x 
and oj = uj r + ij denote a physical parameter, the wave 
number in the x direction, and the angular frequency, 
respectively. 

Linearizing the above MHD equations (p]|4|), one ob- 
tains 



k x V x 8v x - i——Svy + - — Sp + k x h z + 



jSv„ = k v V v Sv„ 



dy 
i d 



2n 



n 



i dB 0x 
n dy 



n 



u5v z = k x V x 6v z - k x ^-b z + — ^^-h y , 
no n dy 

u)5p = k x V x 5p - iSvy--^ + TP k x Sv x - irP — — }L , 
dy 5y 
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Discretizing the spatial derivative of the Fourier am- 
plitude in the y direction with boundary conditions at 
y = ±yt far away from the shear layer, the equations 
cans be cast into the form of a eigen value problem as 



(7) 



from which we can obtain the eigen values of w, whose 
imaginary part is the growth rate, and the corresponding 
eigen functions for any k x . To solve th is eigen value 
problem, we use used the QR algorithm ()Francis![l96lL 
11961 IKublanovskavalH96l . 

4.3. Adopted initial parameters 

Solving the above described problem requires to know 
the initial conditions of the set-up. Most of them can be 
derived from observations, or estimated. The adopted 
physical conditions are summarized in Table [TJ The den- 
sity and temperature (nj, Tj) of the molecular cloud have 
been discussed in BMC. According to their results we 
adopt ni = 10 4 cm' 3 and T/ = 20 K. For the Hll re- 
gion, we consider a typical temperature Tn = 10 4 K, 
and rin — 20 cm -3 . The velocity sheer is tak en to 
be V = 10 km.s- 1 following iRov fc Kunthl (pMh , also 
typical for such environments. Magnetic field strength 
has be en measured in t he Or ion nebula by IBrogan et al.1 
(|2005f) and I Abel et al.1 (|2004f) and found to vary between 
5nT for the Trapezium region and and 25 nT in the 
Veil region. The Ripples are likely situated between 
the Trapezium and the Veil so we adopt a conserva- 
tive value of B = 20 nT. The value of 8 cannot be 
determine independently so we have considered various 
values between and 90° (with symmetrical results in 
the -90 to 0° range). Finally, BMC derived a gravita- 
tional field at the cloud surface g — 3.5 x 10~ n m s~ 2 . 
We can compare the potential to kinet ic energy using 
the R ichardson number Ri = g L/(Vq) 2 (|Chandrasekharl 
[HHI). The value of L can be estimated to be the thick- 
ness of the photodissociation region (PDR) measured by 
BMC, which is essentially the region where the gas is 
converted from fully neutral to fully ionized and where 
the temperature changes from a few 10 K to a few 1000 
K (|Tielens fc Hoellenbachl [l985l) . L is typically 0.01 pc 
(BMC), hence Ri ~ 10 -3 which implies that gravity is 
ineffective and it is no further included in the calcula- 
tions. 

5. RESULTS 

5.1. Linear phase of the KHI 

The results of the linear analysis concern the linear 
growth of the instability. They are presented in Fig. [3] 
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2. Schematic representation of the configuration adopted 
initial conditions of the interface between the H II region 
molecular cloud. 



and Fig. @] Fig. [3] shows the influence of the value of 8, 
the inclination of the magnetic field, on the normalized 
growth rate jL/Vq and on the normalized wavenumbcr 
k x L. The growth rate decreases with increasing angle but 
for —25° < 8 < 25° the growth rate is non-zero, imply- 
ing that the interface is KH unstable. In this range of 
values for 8, the wavenumber k x is almost constant, show- 
ing that it does not depend on magnetic field orienta- 
tion. Fig.|4]shows the evolution of the normalized growth 
rate 7L/V0 as a function of the normalized wavenumbcr 
k x L, for 8 — 0. The most unstable mode corresponds to 
k x L = 0.56, and typically, for 0.1 < k x L < 1 the growth 
rate is high and corresponds to an e-fold timescale of less 
than 10 5 years, short compared to the lifetime of an OB 
association (^10 Myrs). Therefore, it is expected that 
instabilities with 0.1 < k x L < 1 appear in star- forming 
regions. Using L = 0.01 pc, imposes that the wavelength 
of the instability Xkh will range between 0.06 and 0.6 
pc. This number also places limits on the detectability 
of KHIs: at a distance of 1 kpc this is an angular size of 
4-40", and at 10 kpc this is 0.4-4". Hence, KHI struc- 
tures are expected to be of small angular size and can 
only be observed with high angular resolution telescope 
and/or in nearby regions of massive star formation like 
Orion. 

5.2. Saturation of the KHI 

iMatsumoto fc Sekil (|2010f ) studied in details the 2 di- 
mensional evolution of the KHI for conditions similar to 
those presented here and found that the growth of the 
instability leads to saturation. This results in the forma- 
tion of a turbulent layer where the two fluids are mixed, 
over a time-scale of the order of t sat ~ 200L/Vb. Hence, 
using the observed value for L (Table [T]) this results in 
a saturation timescale t sat ~ 2 x 10 5 years. It is impor- 
tant to realize that even if t sat is short there is always 
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Table 1 

Physical parameters 



Parameter 


Symbol 




Reference 


Observational 


Heliocentric distance 


d 


414 pc 


(1) 


Neutral gas density 


m 


10 4 cm" 3 


(2) 


Ionized gas density 


nu 


20 cm -3 


(2) 


Neutral gas Temperature 


Ti 


20 K 


(3) 


Ionized gas Temperature 


T n 


10 4 K 


(4) 


Velocity sheer adopted here 


Vo 


10 km.s" 1 


(5) 


Gravitational field 


9 


3.5 xHT 11 m.s~ 2 


(2) 


Magnetic field strength 


B 


20 nT 


(6) 


Magnetic field orientations 


e 


0-90° 




Instability wavelength 


\ 


0.11 pc 


(2) 


Width of the sheered layer 


L 


0.01 pc 


(2) 


Linear regime length 




0.3 pc 


Fig. 1 


Derived 


values (for 


maximal growth rate) 




Instability growth rate 


7 


2.3xl0- 13 s- 1 




Instability phase velocity 


v+ 


3.6 km.s" 1 




Instability saturation timescale 


tsat 


2xl0 5 yrs 




Size of mixing layer after t sa t 




0.5 pc 




Distance travelled before saturation 


Lsat 


0.74 pc 





(1) From Mcntcn ct al. (2007), (2) From BMC, (3) Typical for molecular clouds (Ticlcns 2005), 
(4) Ty pical for Hn regions (Ticlcns 2Q0D> (5) From lRov & Kunthl (1991), (6) Based on I Abel et al.l 
(2004) 
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Figure 3. Left: Normalized growth rate ■yL/Vo, of the most un- 
stable mode, as a function of the orientation of the magnetic field. 
Right: Normalized spatial wavenumber k x L of the most unstable 
mode as a function of the orientation of the magnetic field. 

a part of the instability that remains linear. Accord- 
ing to the results of our linear analysis we find that the 
KH mode travels along the boundary layer with speed of 
Vcj, = 0.36 x Vo, that is ~ 3.6 km/s. Hence, we can define 
L sat , the spatial scale before the instability has saturated 
by L sat = t sa t and find L sat ~ 0.74 pc. These theoret- 
ical results have several implications. First, the timescale 
for saturation is short compared to the lifetime of an OB 
association, so KHIs will saturate and will be a source 
of turbulence. Secondly, the scale size over which it is 
possible to see the linear regime is at maximum 0.74 pc. 

6. THE RIPPLES AS A KHI 

Based on the theoretical results described above, we 
investigate in more details if the Ripples can be inter- 
preted as an occurrence of the linear phase of the KHI. 
The observed value of A for the Ripples is 0.1 pc, which 
fall in the range defined from the theoretical investigation 
(0.06 < Xkh < 0.6 pc). The Ripples seem to preserve a 
very periodic structure, suggesting linear regime, over a 
distance Lu n — 3A = 0.3 pc (Fig. 1). The following 2 bil- 



0.02 t 



0.01 




Figure 4. Normalized growth rate jL/Vo versus normalized spa- 
tial wavenumber k x L of the instability. The right axis translates 
the value of the growth rate into an e-fold timescale for Orion. 



lows, instead, start showing some chaotic structure sug- 
gesting the beginning of the saturation of the instability. 
After about 5A, the periodic structure has disappeared. 
Hence, the travelled distance in linear regime Lu n for 
the Ripples is, consistently with the KHI model, smaller 
than the maximal theoretical value of L sat — 0.74 pc dis- 
cussed in Sect. [SJ Therefore, we argue that the observed 
evolution of the Ripples structure results from the mo- 
tion of the KH wave towards saturation, from left to 
right in Fig. 1. Our last remark concerns the impor- 
tance of the magnetic field. As mentioned above the 
Ripples can only result from a KHI if they are in a re- 
gion of small 8. Orientation of the magnetic field lines i n 
Orion have been measured by e.g. iHoude et al.l (|2004l) ; 
iPoidevin. Batien fc Matthews! ({2010) using polarimetry. 
Unfortunately, this does not cover the Ripples region 
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and, in addition, performing such measurements at the 
arc-second scale remains challenging. We can only stress 
the importance -often neglected- that magnetic fields 
have in shaping the interstellar medium, in this case be- 
cause they can stabilize interfaces between H n regions 
and molecular clouds. 

Altogether, we conclude that the study presented in 
this letter brings additional evidence that the Ripples 
result from a KHI. This raises one question however, 
which is "why do we only see one occurrence of the KHI 
in Orion?" This is perhaps because only in this region 
are the conditions (e.g. magnetic field, velocity flow) 
favorable at the moment we observe Orion. Over the 
lifetime of the region however, this may occur a high 
number of times. Ripples have indeed been observed 
rec ently in an o ther s tar forming region (Cygnus OB2) 
by iSahai et al.l (|2012l ) , suggesting that the structure in 
Orion is not unique. It is also possible that more of these 
structures exist in Orion, but given their small angular 
size or unfavorable orientation on the plane of the sky 
they remain undetectable (see Sect. [5]). 

7. THE KHI AND CHEMICAL MIXING IN STAR-FORMING 
REGIONS 

7.1. Saturation of the instability: turbulence and 
chemical mixing 

We have discussed in Sect. [5] the saturation of the 
KHI towards a turbulent regime. This may have an im- 
portant role in chem ical mixing in star forming regions. 
IRov fe Kunthl (fl99lh were the first to recognize the im- 
portance of the KHI in chemical mixing of the interstellar 
medium. They studied the influence of KHIs by defining 
the e-fold timescale assuming a spatial wavelength of 100 
pc (the size of a large Hn region), and from this derived 
a timescale of 1.5 x 10 6 years, which is smaller than the 
lifetime of an OB association. However, their model was 
based on an ideal hydrodynamical case with no magnetic 
field, for which the wavelength of maximal growth rate 
cannot be determined. This is partly due to the fact 
that no observational evidence of KHIs existed at the 
time. The results presented in this Letter using a more 
detailed model and guided by direct observations, are 
clearly incompatible with a maximal gr owth rate cor- 
respon ding to A = 100 pc. In addition, IRov fe Kunthl 
(1991) used the e-fold timescale as a measure of mixing 
timescale, which is not appropriate. Instead, we can use 
the re sults found here and those of [Matsumoto fe Sekil 
(2010) to obtain some general appreciation of the effi- 
ciency of the KHI mixing in star-forming regions. First, 
as mentioned above, full mixing of the fluids occurs after 
t sat which we have found to be of the order of a few 10 5 
years. Again, this is short compared to the lifetime of 
an OB association so the process will be efficient. After 
this time, the size of the mixed layer is L m i X = 50 x L 
(jMatsumoto fe Sekil 120101 Fig 91. that is 5 pc. All 
in all, in agreement with lllov fe Kunthl (|1991l ) (although 
making different hypotheses), we conclude that the KHI 
can be an efficient mechanism to mix chemical elements 
in the interstellar medium. 

7.2. Further implications: 26 Al in the solar system 

It is believed that low-mass star formation is triggered 
in the over-dense shell of molecular cloud that lie around 



Hn re gions (see e.g. [Deharveng et al.ll2010f) . In a recent 



Gounelle fe Mevnetl (|2012ft argue that the Solar 



paper, 

system may have formed in such an environment, based 
on the abundances of short lived radionuclides found in 
meteorites. In particular, they propose that 26 Al was 
brought to the forming solar system by the w inds of 
a massive star (see also iMontmerle et al.l |2007|) rather 
than by Supernovae. This requires that a nearby mas- 
sive star (M* > 32Mq) injected 26 Al during a few million 
years, and that this element was then well mixed with the 
Hll gas, and eventually that the Hn gas was mixed with 
the surrounding molecular shell efficiently. In Orion at 
least, the gas from the wind seems to be well m ixed with 
the Hn region as shown bv lGudel et~aT1 (|2008l ). The ef- 
ficiency of H II mixing with the molecu lar shell was not 
evaluated bv I Gounelle fc Mevnetl <|2012f ) . However, they 
derive the time during which the molecular shell has 
to be enriched before the solar system starts to form. 
This value ranges between 0.65 and 6.2 Myr. The mix- 
ing timescale we have derived here (t sat = 2 x 10 5 yrs) is 
smaller than t+ so that indeed mixing by the KHI is effi- 
cient to enrich the molecular shell with 26 Al. Hence, the 
KHI (and possibly other instabilities) could have played 
an important role in the transport of 26 Al to the forming 
Solar system. 

8. CONCLUSION 

We have shown that the KHI develops rapidly at the 
H II molecular cloud interface in conditions like Orion 
(which are representative of many massive star forming 
regions) . After travelling at the surface of the cloud dur- 
ing a time of a few 10 5 years and over a maximum dis- 
tance of ~ 0.74 pc, the instability reaches saturation. 
Hence, as suspected, the KHI is probably a significant 
mechanism to generate small scale (< 1 pc) turbulence 
in molecular clouds near massive stars. In addition, since 
the H II region is contaminated by the chemical elements 
injected by massive stars winds, the KHI may be a rele- 
vant process to bring these elements inside the molecular 
cloud, in regions where planetary systems around young 
stars are formed. This could have been at play to trans- 
port 26 Al to the nascent solar system. Periodic structures 
corresponding to the linear phase of the KHI (like the 
Ripples) should be relatively widespread in star-forming 
regions, for instance on the surface of m olecular globules 
as reported recently (jSahai et al.ll2012|) . However, these 
structures are expected to be small (few arc-seconds) and 
hence require high angular resolution observations to be 
identified. 
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